Decay of the Mixed States 
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We study the classical escape from local minima for 2d multi-well Hamiltonian systems, realizing 
the mixed state. We show that escape from such local minima has a diversity of principally new 
features, representing an interesting topic for conceptual understanding of chaotic dynamics and 
applications. 
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The escape of trajectories (particles) from localized re- 
gions of phase or configuration space has been an im- 
portant topic in dynamics, because it describes the de- 
cay phenomena of metastable states in many branches 
of physics: chemical and nuclear reactions, atomic ion- 
ization, nuclear fusion and so on. This problem has the 
rich history. Almost a century ago, Sabine [l[ consid- 
ered the decay of sound in concert halls. Legrand and 
Sornette Q have shown that this problem is equivalent 
to the escape one: a small opening of width A for es- 
cape must be identified with J a(S)ds, where a(s) is the 
absorption coefficient at position S of the container (bil- 
liard) boundary, a(s) = 1 over the width of window and 
a(s) = elsewhere. Szepfalusy and Tel @ connected 
escape problem with problem of chaotic scattering. 

Exponential decay is a common property expected in 
strongly chaotic classical systems @, H, Q. Let us con- 
sider as an example ; 4] point particles bouncing elasti- 
cally off the walls in a rectangular box. The system is 
allowed to decay by providing a small window in one of 
the box walls through which particles can escape. As is 
well known, motion of particles in a rectangular billiard 
is regular: two independent integrals of motion are abso- 
lute values of momentum projection on the billiard walls. 
The trajectories of particles become chaotic if a circular 
scattering center is placed somewhere inside the box. 

For the chaotic case simple consideration leads to the 
exponential decay. The number of particles leaving per 
time interval is given by 



Ap{t) / d 2 P pe„ = -2Ap(t)p 2 S P (1) 



dN 
~dt 

Here p is absolute value of the momentum, e„ is a unit 
vector normal to the opening in the surface, and inte- 
gration in momentum space is taken over a circular ring 
with radius p and infinitesimal width Sp. Function p(t) 
is the phase space density, which for ergodic motion is 
only a function of time. In our case 



P(t) = 



2irp6pA c 



(2) 



where A c is total coordinate space area available. Insert- 
ing © into (dJ yields 



N(t) = N(0)e~ 



pA 

7tX 



(3) 



Analytically calculated decay constant a is in a good 
agreement with the graphically extracted value. 

Exponential law at extremely long times turns into the 
power law typical for decay of regular systems. One pos- 
sible mechanism for generation of power tails is the ef- 
fect of "sticking" of the chaotic orbits to outer bound- 
aries of stability islands Q, or a very similar effect, con- 
nected with the existence of marginally stable periodic 
or "bouncing ball" orbits. Although some qualitative 
models, which show how the algebraic tail emerges, were 
introduced in Q, no critical conditions for the distinct 
decay laws were formulated in terms of the billiard geo- 
metrical constrains. Experimental escape of cold atoms 
from a laser trap of billiard type with a hole was studied 
in0,[l|. 

Transition from the billiards to potential systems sub- 
stantially broadens number of possible applications of the 
escape problem, but from another hand significantly com- 
plicates the problem. Of course, the one-well case is the 
simplest one. Zhao and Du reported a study on the 
escape rates near threshold of Hcnon-Hciles potential 
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Simulations performed by the authors show that the es- 
cape of Henon-Heiles system at energy, slightly exceed- 
ing the saddle one, follows exponential law similar to 
the chaotic billiard systems. They derived an analytic 
formula for the escape rate as function of energy. The 
derivation is based on the fact that the phase space of the 
considered potential (as well as for the billiards) is prac- 
tically homogeneous near the saddle points. It should be 
noted that in such case all trajectories with energy higher 
than the saddle one leave the potential well in finite time. 
The only problem to solve is to determine the probability 
of particle to escape from the well in unit time interval. 

In contrast to billiards, generic potential systems have 
essentially inhomogeneous phase space structure. We in- 



FIG. 1: Level lines for potential D 5 © (a) and QO @ (b). 
Gray color shows initial distribution of particles for E = 2Es- 

tend to study the particles escape from the local min- 
ima in the case when the phase space contains macro- 
scopically significant components of regular as well as of 
chaotic type. Such possibility is realized in multi-well 
potentials. 

The principal peculiarity of the regularity-chaos tran- 
sition in multi-well potentials lies in the existence of 
different critical energies for different local minima, ft 
means that in such potentials at one and the same energy 
in different local minima may exist different dynamical 
regimes (either regular or chaotic). Such kind of dynam- 
ics in multi-well potentials, when at some energy the ratio 
of chaotic trajectories in certain local minimum signifi- 
cantly differs from that ratio in other minima, is called 
the mixed state [l2T |. 

We demonstrate the mixed states on two representative 
examples: the lower umbillic catastrophe D5 potential 

m 

U D5 (x,y)=y 2 (x + 2a)+(j-l) (5) 

for a = 1.1 (fig[T]a) and the potential of quadrupole os- 
cillations (QO) of atomic nuclei [HI 

U Q o(x, y) = X ^±f- + V ~ y + (x 2 + V 2 ) 2 (6) 

for W = 18 (figfflb). The potential D 5 (JSJ has only 
two local minima and three saddles and it is the simplest 
potential, where the mixed state is observed. 

FigH] presents the Poincare sections for different en- 
ergies, demonstrating evolution of dynamics in different 
local minima. At low energies motion has well-marked 
quasiperiodic character for both minima (fig|2la,e). As 
energy grows, gradual regularity-to-chaos transition is 
observed. However changes in features of the trajecto- 
ries, localized in certain minima, are sharply distinct. For 
the left minimum, already at about half saddle energy, 
significant fraction of the trajectories becomes chaotic 
(fig[2jb,f), and at saddle energy practically all initial con- 
ditions produce chaotic trajectories (figH]c,g). In right 
minimum under the same conditions motion remains 
quasiperiodic up to the saddle energy (further we will 



FIG. 2: Poncare surfaces of section y — in the (x,p x ) plane 
for the potentials D5 ([5| (a-d) and QO (J6]) (e-h) for different 
energies (from left to right): E < Eg, E > E s , E = Es, 
E — 2Es- Solid line delimits the classically allowed region of 
phase space. 

call it "regular local minimum" for simplicity). More- 
over, at energies significantly higher than the saddle en- 
ergy (see fig HI , the phase space structure preserves divi- 
sion on chaotic and regular components (figHJd,h). The 
latter is localized in the part of the configuration space 
which corresponds to regular motion at energies below 
the saddle. Earlier we have shown that the mixed state 
opens new possibilities for investigations of quantum sig- 
natures of classical stochasticity [TH, Ht| . Aim of the 
present work is to study the classical escape from sepa- 
rated local minima, realizing the mixed state. We show, 
that escape from such local minima has all the above 
mentioned properties of decay of chaotic systems, and 
also a diversity of principally new features, represent- 
ing an interesting topic for conceptual understanding of 
chaotic dynamics, and for the applications as well. We 
are interested only in the "first passage" effects, leaving 
aside the problem of dynamical equilibrium setup for the 
finite motion (for example, in QO potential). In is impor- 
tant to stress that though we study the process of escape 
from a concrete local minimum, the over-barrier in the 
case of mixed state has a specific memory: general phase 
space structure at super-saddle energies is determined by 
the characteristics of motion in all other local minima. 

We carried out numerical simulation and analytical es- 
timates of trajectories escape in the potentials D5 and 
QO through the hole over the saddle point. Results of 
the escape problem for systems with multi-component 
phase volume (regular and chaotic components) essen- 
tially depend on choice of ensembles of initial conditions 
for dynamical variables. Fig|3] presents the normalized 
particle number N(t)/N(t = 0) for 10 6 initial conditions, 
uniformly distributed inside the right minimum in the 
potential D§ (x > 0) and peripheral minimum in the 
QO potential (x > 1/12) together with the typical tra- 
jectories and Poincare sections. The results for different 
potentials are evidently similar and have such character- 
istic features: 

• At times t — > 00 decay law saturates at 
N(t -> 00) = p^N 
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FIG. 3: Decay law for mixed states in the D5 (a) and QO (b) 
potentials. Solid lines — numerical simulation for E/Es = 
1.1, 1.5, 2.0; dotted and dashed lines — theoretically obtained 
exponential and linear decay laws respectively, zoomed on 
the inset figure in the lower left corners. Other inset figures 
represent the typical trajectories and Poincare sections for the 
three different types of initial conditions: linearly escaping, 
exponentially escaping and non-escaping. Circles show the 
joining points between the linear and the exponential decay 
laws at critical time t = r. 



where p( ne ' is equal to relative phase volume of 
"never-escaping" trajectories, which represent the 
regular trajectories, completely localized inside the 
considered minimum. All such trajectories, there- 
fore, have infinite escaping times. 

• For t > t(E) the decay law has exponential form 

N(t)/N = P {ne) + pto e - aM Q- T ) (7) 

where p^ represents relative number of exponen- 
tially escaping particles. 

• For t < t(E) the decay law is linear 

N(t)/N Q = 1 - a®t. (8) 

We should stress, that ((8]) is in no way a linear approxi- 
mation of (O for small t: in general p( ne ' + pW e~ a{ " ) r 7^ 1 

and a( e )p( e )e~ Q<e>T 7^ . Instead, from the condition 
of smooth joining of curves ([8} and (|7|) in the transition 
point t = t we obtain 



,W = 1 



where p"> = aS 'T is the relative number of linearly escap- 
ing particles. Moreover, already on time scales t < t(E) 
linear decay law ([8|) is apparently different from its ex- 
ponential analogue 



(ne)N (-a«'t) 



N(t)/N{t = 0) = p (ne > + (1 - p < - ne >)e 

(see the inset on figE]). 

As one can see from the inset Poincare section on fig[3] 
both the chaotic and regular trajectories contribute to 
linear escaping regime ([8j). because for sufficiently small 
times t < t chaotic and regular motions are not yet dis- 
tinguishable. Up to transient time t = r all quasi-one- 
dimensional regular trajectories, oriented along the re- 
axis, already escape and for t > r the escape of remaining 
chaotic particles follows exponential law ([7|). The parti- 
cles escaping the last show already mentioned sticking 
phenomenon (see the inset Poincare sections on fig [3]). 

The transient time t(E) in fact coincides with the pas- 
sage time of the longest one-dimensional path from the 
opening to the opposite wall of the potential well and 
back (see figE]). For the potentials D5 and QO corre- 
sponding theoretical estimates read (we assumed m = 1) 



td 5 (E)=2 f 



I 



TQO (E) = 12 ) k K 



m 1^ 




where K(k) is the complete elliptic integral of the first 
kind and E$ = 1/12 4 is the saddle energy in the QO 
potential for W = 18 (for the D 5 potential Es = 1). 

Theoretical estimates for the escape rate were obtained 
by averaging the escape probability over the opening [Til ] : 



(E) = p{E) dy d9\p\ cos 8 



where xs is the coordinate of the saddle point and p(E) 
is the normalized particle density: 

P(E) 



2irA(E) ' 



where A(E) denotes area of the classically allowed region 
inside the well: 



ME) 



dxdyQ (E-U(x,y)) 



x>xs 



p (ne) _ p (l) 



Such p(E) corresponds to uniform distribution of initial 
conditions on the energy surface H(p,q) = E. 

For the potentials D$ and QO the explicit formulae are 
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x + 2a 



Aqo(E) = i 



«VK + 4) 2 - 7 x 



- (e 2 - 1)' 



[(£ + 4) 2 -7p 



where £ = \/Es{x — xg). 

Finally, the general expression for the escape rate is 



a(E) 



(P> 



ttA(£;) ■ 

For our case (p) = J d!y|p|. In the case of billiards with 

small opening p = const =>• (p) = pA, and we recover 
expression (J3|) . For the potentials D 5 and QO we get the 
results in closed form: 



uqo(E) = 



E - 1 



2y/EA Ds {E) 



12ttA qo (E) 



2E 



where s = E — E$ + 1/256, K(k) and E(k) are the com- 
plete elliptic integral of the first and second kind respec- 
tively. 

In order to obtain (E) we correct A(E) subtracting 
the relative phase space occupied by the non-escaping 
particles 



A (l \E) =A(E)(1 -p(" e )) 



,(0 



(E) 



a{E) 
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FigE] demonstrates good agreement between our theoret- 
ical and numerical results for wide energy range. 

Fraction of the non-escaping particles p( ne ) coincides 
with the relative phase space volume of trajectories, lo- 
calized in the regular minimum, which may be well es- 
timated by the relative area of the stability island p^ 
on the Poincare section (fig©. Calculation of relative 
area of regular island in the Poincare section was per- 
formed by the following scheme. First, the island bound- 
ary was determined by numerical integration of equations 
of motion and then the interior area was calculated. Fur- 
ther the obtained area was divided on the entire area 



E/Es 



FIG. 4: Correlations between relative area of the stability 
island p {si \E) and non-escaping particles ratio p'" e ': empty 



squares - p^{E), filled squares - p [ ^ e> (E), empty circles 



(ne) , 



d${E), filled circles -p%g(E) 



(ne) , 



of classically allowed motion, defined by the conditions 
x > and p 2 > 0. While the phase volume itself is 
4-dimensional, the stability island in Poincare section is 
2-dimcnsional, so we cannot expect absolute coincidence 
of the corresponding measures. However, the calcula- 
tions show very close correspondence between them (see 
figlll)- Therefore, numerical analysis of Poincare sections 
together with our theoretical results gives all information 
necessary to predict the escape dynamics in independent 
way. 

In summary, we have considered classical escape from 
separated local minima in two representative 2D multi- 
well potentials, realizing the mixed state. We have found 
that escape from regular minima contains a number of 
new features. The most important among them are the 
following: 

1. Decay law saturates at long time ranges. 

2. On small time scales there exists a linear segment, 
which is not connected with linear approximation 
to the exponential decay law, observed in chaotic 
systems with homogeneous phase space. 

3. Fraction of particles, remaining in the well, is deter- 
mined by relative phase volume of the regular com- 
ponent, which in its turn monotonically decreases 
with growth of energy. 

It was shown that the linear segment of the decay law 
is generated by the quasi-one-dimensional trajectories, 
oriented perpendicular to the opening, and the transient 
time of the linear-to-exponential regime lies in perfect 
agreement with the analytical estimates. 

We should note that we devote main attention to es- 
cape from the regular local minima because the specifics 
of the mixed state manifests only in them. However let 
us remind that in the case of mixed state the phase space 
structure at super-saddle energies is determined by dy- 
namical characteristics in different local minima of whole 
potential energy surface. 
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Above mentioned peculiarities of the escape problem 
may found practical application for extraction of required 
particle number from atomic traps. Changing energy of 
particles trapped inside the regular minimum, we can 
extract from the trap any required number of particles. 
Problem of particle energy changing in the potential well 
may be solved by introduction of small dissipation. Ob- 



tained results may present an interest also for description 
of induced nuclear fission in the case of double-humped 
fission barrier. Revealed peculiarities must manifest also 
in over-barrier dynamics of wave packets, initially local- 
ized in the regular minima. 
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